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ABSTRACT 


The shape of the spectral energy distribution of active galaxies in the EUV- 
soft X-ray band (13.6 eV to 1 keV) is uncertain because obscuration by dust 
and gas can hamper our view of the continuum. To investigate the shape of the 
spectral energy distribution in this energy band, we have generated a set of pho- 
toionization models which reproduce the small dispersion found in correlations 
between high-ionization mid-infrared emission lines in a s^ple of hard X-ray 
selected AGN. Our calculations show that a broken power-law continuum model 
is sufficient to reproduce the [Ne V]i 4 . 32 ;ini/[Ne III], [Ne V]24.32/im/[0 IV]25.89/im 
and [O IV] 25 . 89 /im/[Ne III] ratios, and does not require the addition of a “big 
bump” EUV model component. We constrain the EUV-soft X-ray slope, otj, to 
be between 1.5 - 2.0 and derive a best fit of ~ 1.9 for Seyfert 1 galaxies, 
consistent with previous studies of intermediate redshift quasars. If we assume a 
blue bump model, most sources in our sample have derived temperatures between 
Tbb = 10®’^® K to 10®'^ K, suggesting that the peak of this component spans a 
large range of energies extending from ~ A600 A to A1900 A. In this case, the 
best fitting peak energy that matches the mid-infrared line ratios of Seyfert 1 
galaxies occurs between ~ A700-A1000 A. Despite the feict that our results do 
not rule out the presence of an EUV bump, we conclude that our power-law 
model produces enough photons with energies > 4 Ry to generate the observed 
amoimt of mid-infrared emission in our sample of BAT AGN. 


Subject headings: AGN: general - galaxies: Seyfert - X-rays - IR 
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1. Introduction 

Active galactic nuclei (AGN) produce an extraordinary amount of energy which is 
thought to be primarily thermal emission from an optically thick accretion disk fueling a 
supermassive black hole (e.g., Peterson et al. 2004). Accretion can produce a characteristic 
feature in the ultraviolet continuum of most AGN known as the “big blue bump” (bbb, 
Shields 1978; Malkan & Sargent 1982) and it has been argued that its energy peak lies 
around the Lyman limit (e.g., Zheng et al. 1997; Telfer et al. 2002; Shang et al. 2005). The 
problem is that in AGN the soft X-ray-extreme-ultraviolet spectral energy distribution 
(SED) can be hidden because of dust and gas, thus making a direct view of the full energy 
distribution impossible. A common approach is to study high redshift quasars where the 
rest-frame extreme-ultraviolet (EUV, A ~ 100-912 A) can be observed from the ground. 
However, for the local universe at z < 0.05 we need to rely on space-borne telescopes and/or 
other indirect methods to uncover the shape of the SED in this particular range of energies. 

Within the paradigm of the unified model of AGN (e.g., Antonucci 1993) both Type 1 
and Type II AGN are intrinsically the same, with their observed difference being due to 
the visibility of the central engine. An optically-thick toroidal structure surrounds the 
central source, located at distances just beyond the broad line region (BLR). Emission 
that arises further out in the narrow line region (NLR), hundreds of parsecs in extent, is 
expected to be fairly unobstructed and isotropic, i.e., independent of our viewing angle. 
However, absorption, dust extinction and star formation may hamper the use of soft X-ray, 
optical and infrared observables on these larger- than-torus scales (e.g., Mulchaey et al. 
1994; Kraemer et al. 2010). 

One way to bypass the above limitations is to focus on high-ionization (E > 40 eV), 
mid-infrared emission lines. To be produced, they require photons with energies between 
~40-100eV (EUV-Soft X-ray), and so we can assume that they are primarily associated 
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with the AGN, in other words, they axe less contaminated by star formation processes. 
Also, they are less sensitive to dust extinction than optical lines. In this regard, it is 
possible to fit these observed lines via photoionization models in order to constrain the AGN 
continuum in this energy band (e.g., Alexander et al. 1999, 2000). Among the mid-infrared 
properties of Seyfert galaxies, the correlations between high-ionization emission lines, 

[Ne III] A15.56 /im, [O IV] A25.89 /xm and [Ne V] AA14.32/24.32 /Ltm, are strikingly tight, 
with a dispersion of less than ~ 0.2 dex, both in flux and luminosity (see Gorjian et al. 
2007; Dudik et al. 2009; Goulding & Alexander 2009; Melendez et al. 2008b; Weaver et al. 
2010; Pereira-Santaella et al. 2010). The [Ne V] emission, given the high energy require to 
be produce (E >98.9 eV), has been unambiguously identified with the presence of an AGN 
(see Satyapal et al. 2007; Abel Sz Satyapal 2008; Bernard-Salas et al. 2009; Weaver et al. 
2010), therefore, the strong correlation between this emission and the lines produced from 
transitions of [Ne III] (E > 41.0 eV) and [0 IV] (E > 54.9 eV) suggest that these lines 
are produced primarily by the AGN. Moreover, Melendez et al. (2008a) probed the utility 
of high-ionization mid-infrared lines as an indicator of AGN power. They found a tight 
correlation, with a dispersion of ~ 0.3 dex, in Seyfert 1 galaxies, between the [O IV] and 
the X-ray 14-195 keV continuum luminosities. A weaker correlation was found for Seyfert 
2 galaxies, which we believe is due to the effect of Compton Scattering in the 14-195 keV 
band in some Seyfert 2 galaxies. This result has since been confirmed in similar studies 
(e.g., Rigby et al. 2009; Diamond-Stanic et al. 2009; Weaver et al. 2010). 

Our aim is to recover the hidden portion of the spectral energy distribution of AGN 
between the EUV and the soft X-ray (13.6 eV to 1 keV). To do this we have performed 
extensive photoionization calculations to characterize the small dispersion foimd in the 
correlations between high-ionization, mid-infrared emission-line ratios. We compare our 
calculations with emission line ratios from Spitzer/IRS spectra of the Swift Burst Alert 
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Telescope (BAT) AGN sample^ of the local universe, presented in Weaver et al. (2010). All 
of our calculations have been carried-out by the photoionization code CLOUDY, version 
08.00, last described by Ferland et al. (1998). 


2. The Models 

We have generated two sets of photoionization models to study the intrinsic SED of 
the ionizing continuum. The first assumes a broken power law, as used by Melendez et al. 
(2008a) and sinoilar to that suggested for NGC 5548 and NGC 4151 (Kraemer et al. 1998, 
2000), of the form Fi, oc with a= 0.5 below 13.6 eV, cq from 13.6 eV to 1 keV and 0.8 
at higher energies. To examine specifically the EUV/soft X-ray band, we varied ai from 1.0 
to 2.5. Our second model is a parametrization of the AGN continuum that combines an 
X-ray power law and a “blue bump” component, in the form; 


fv = 1 ^““” exp {-hu/hTBs) exp {-kTin/hv) -f- ai/**, (1) 


where cruv is the low-energy slope of the Big Blue Bump continuum, Tbb is the UV 
bump cut-off temperature and kTin is the infrared cut-off temperature of the big bump, 
corresponding to 9.1 /im. For the X-ray power-law, cxx is the slope of the X-ray component 
and the coefficient a is adjusted to produce the correct UV to X-ray slope (Tananbaum et al. 
1979), a<,a:, as defined by 


/.(2keV) 

/^(2500A) 


403.3““. 


( 2 ) 


^The ^itn/t/BAT covers the whole sky at (l-3)xl0~^^ergs cm“^ s“" and represents a 
complete sample of active galaxies including Compton-thin AGN that were missed from 
previous X-ray surveys in the 2-10 keV band because of their high column densities {Nh ~ 
10^^ cm“^ ). 
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For the “big bump” model we adopted a value for the UV to X-ray spectral slope of 
aox = —1-4 (Zamorani et al. 1981), in good agreement with recent studies of Type I 
AGN (Lusso et al. 2010). We fixed the low-energy slope of the big blue bump continuum 
to = —0.5 (e.g., Francis et al. 1993). For the X-ray power-law, we adopted a value 
of 0 !a, = — 1. This parametrization is similar to that of Mathews & Ferland (1987) and 
Korista et al. (1997). Finally, for this model we considered five values for the UV biunp 
cutoff temperature, from Tbb = 10®'^® K to 10® ° K. Figure 1 compares the SED for our 
different photoionization models. 


To investigate the physical conditions in the emission line regions for the [O IV], [Ne III] 
and [Ne V] lines, we started with a simple, single-zone, constant density model. The logs of 
abimdances of Oxygen and Neon relative to H by number are -3.31 and -4.0, respectively, 
where the O abundance is from Allende Prieto et al. (2001) and the Ne abundance is from 
Holweger (2001). We considered hydrogen column densities of lO^^cm"^ and lO^^cm"^. 
Despite the fact that dust is likely mixed in with the emission-line gas in the NLRs of 
Seyfert galaxies (e.g., Kraemer Sz Harrington 1986; Netzer k Laor 1993), the tightness of 
the mid-infrared correlation suggest that extinction in the mid-infrared, on average, is not 
an important bias in our sample of X-ray selected AGN (Weaver et al. 2010). 


We generated a grid of photoionization models varjdng the total hydrogen number 
density, (n^), and the ionization parameter, U, where the ionization parameter U is defined 
as (see Osterbrock k Ferland 2006): 


„ 1 r n , <3(g) 

AwF^criff hv iirR^cnn ’ 


( 3 ) 


where R is the distance to the cloud, c is the speed of light and Q{H) is the flux of ionizing 
photons. In the next section we will Compare these models with the line ratios observed 
from the BAT AGN sample. 
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3. Comparison with Observations 
3.1. Density in the Emission Line Regions 

The use of emission lines that are close in ionization potential but different critical densi- 
ties can be used as a density diagnostic of photoionized plasmas (e.g., Osterbrock &: Ferland 
2006). In particular, the mid-infrared [Ne V] ratio, [Ne V]i 4 . 32 /im/[Ne V] 24 , 32 /im» has been 
used as a density diagnostic in active galaxies (e.g., Dudik et al. 2007; Tommasin et al. 2010) 
with the caveat that a wavelength-dependent mid-inftared extinction can affect this ratio, 
as suggested by the number of Seyfert galaxies that have [Ne V] ratios below the theoretical 
low density limit (e.g., Tommasin et al. 2010; Weaver et al. 2010; Pereira-SantaeUa et al. 
2010). However, the mid- infrared extinction curve derived by Weingartner & Draine (2001) 
and Li & Draine (2001) (e.g, see Figure 16 in Li & Draine 2001) fails to predict the required 
dust extinction at Ai 4 . 32 /im/NH and A 24 . 32 /im/NH to explain the observed [Ne V] ratios 
that fall below the low density limit in the BAT sample (Weaver et al. 2010). Therefore, 
since the relation between the [Ne V] ratio and extinction is not fully understood, we 
minimized the possible effect of reddening toward the NLR by focusing on the [Ne V] 
ratio observed in Seyfert 1 galaxies as detailed multi- wavelength studies on Se 3 dert galaxies 
suggest more extinction of emission from the NLR in Se)dert 2 galaxies than in Seyfert Is 
(e.g., Melendez et al. 2008a; Kraemer et al. 2011). Following the work on the mid-infrared 
properties of the BAT sample of AGN (Weaver et al. 2010) we found an average neon ratio 
observed in Seyfert 1 galaxies of [Ne V]i 4 . 32 //m/[Ne V] 24 . 32 /im = 0.95 ± 0.28 (see Table 1). 

We used the optimize^ command in CLOUDY to vary and U to reproduce the 
observed neon ratio. From this, and by using our “power-law” model with a.: = 1.5 for a 

^For afull description on the optimization methods see CLOUDY documentation. Hazy 1, 
at http://www.nublado.org/ 
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total hydrogen column density of uh = we obtained a value of Uh = 10^ ®® (cm“^) 

and log 17 = —1.37. For our “blue bump” model, assuming a cutoff temperature of 
Tbb = 10®'^ K such that the energy of the UV bump peaked at ~ 22 eV (e.g., Korista et al. 
1997), we obtained a value of rih = 10^'^^ (cm"®) and log 17 = —1.51. For a toted hydrogen 
column density of Nh = 10®® cm"® we obtained a value of = 10® ®® (cm"®) and 
log!/ = —1.31 for oiu power-law models and, from our “blue bump” model, assuming 
a cutoff temperature of Tbb = 10®'^ K, we obtained a value of Uh = 10® ®^ (cm"®) and 
log 17 = —1.30. In Figure 2 we present the neon ratio versus the total hydrogen number 
density for these models for a range of ionization parameters and EUV/Soft X-ray slopes. 
Figure 2 shows that, close to the observed values, the neon ratio is not strongly sensitive to 
the SED of the ionizing radiation. 

From the different models and continuum approximations we obtained a total hydrogen 
number density of n* « 10® (cm"®). One must note that this value represents the density 
only of the [Ne V] emission-line gas. However, this density is in agreement with the 
values found from the photoionization modeling of the [O IV] emitting region presented 
in Melendez et al. (2008a), suggesting that the [Ne V] and [O IV] emitting regions are 
essentially co-located. On the other hand, this density is smaller than that found by 
Baskin & Laor (2005) for the [O III] emitting region in their single-zone approximation, i.e., 
He ~ 10®-®® (cm"®). However, they found that the observed [0 III] A5007, [O III] A4363 and 
H,? lines can also be fitted using two-zone model comprised of a a dense (ng ~ 10® cm"®) 
inner zone, dominated by the [O III] A4363 emission, and a low density, more extended zone 
(ug ~ 10® cm"®), where most of the [0 III] A5007 arises. Our density value is similar to that 
used in their low-density zone, although the densities in their two-zone approximation are 
based on fixed assumed values, whereas our density is derived from the fit to the observed 
[Ne V] ratio. 
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3.2. The EUV— Soft X-ray slope 

Following the results from the previous section, and in order to keep the number of free 
parameters to a minimum, we used the observed [Ne V] ratios to fix the hydrogen density 
in our photoionization models to n/, = 10^ (cm“®). Next, in order to maximize each of the 
mid-infrared emission lines, we proceeded to find the range in ionization parameter where 
the predicted ionic column densities for the Ne III, Ne V and O IV peak (see Figure 3). 
Results show that the ionization parameter falls in the range —2.0 < logt/ < —1.5 for 
Nh = lO^^cm"^, while it is in a higher ionization range, logf/ ~ —1.0, for Nh = 10^cm“^. 
One must note that the peak of the ionic column density depends on the SED (see Figure 5). 
Also, because we are interested in the shape of the ionizing continuum, we expect to achieve 
the clearest constraints from radiation-bounded models The set of radiation- bounded 
models was obtained by initially assuming a hydrogen column density of Nh = lO^^cm"^. In 
Figure 4 we compared the total Hydrogen density as a function of ionization parameter and 
UV slope for our radiation bounded, “power-law” models. It is clear that a column density 
of Nh = lO^cm"^ is enough to generate an optically thick, radiation-bounded cloud for the 
set of ionization parameters and hydrogen density derived from the previous considerations, 
except for ai = 1.0 where we had to extend the column density to Nh = lO^^cm"^. Using 
this set of ionization parameters, and n/i = 10^ cm“®, we compared our model predictions 
with the observed [Ne V]i 4 . 32 ^m/[Ne III], [Ne V]24.32/im/[0 IV] and [0 IV]/[Ne III] ratios 
from the BAT sample. Overall, these results suggest that the [O IV] and [Ne V] come from 
a higher ionization state and lower density regions than that suggested by the optical [O III] 
lines (Baskin & Laor 2005; Melendez et al. 2008a). 

^Specifically, in a radiation-boimded model, all of the incident radiation has been absorbed 
inside the cloud, i.e., the cloud is optically thick to the ionizing radiation. 
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3.3, The Power-Law model 

As we mentioned before, high-ionization mid-infrared ratios can be used to investigate 
different excitation mechanisms, thus different shapes of the ionizing continuum. They 
provide information about the relative flux density of ionizing photons with energies 
between ~40-100 eV. In Figure 6 we compared the observed [Ne V]i4.32/im/[Ne III] and 
[Ne V]24.32/ira/[0 IV] ratios found in the BAT sample (see Table 1) against those predicted 
by our photoionization calculations. Because of the proximity of wavelengths and equal 
aperture extraction for the emission lines in each ratin'*, this comparison is unaffected 
by dust extinction and aperture effects. Figure 6 shows two sets of “power-law” models, 
corresponding to Njj = 10^* cm~^ and Nh = 10^^ cm~^, for a range of the EUV/soft X-ray 
from Qj = 1.5 — 2.0. We show that, for a column density of Nff = 10^^ cm“^, the cloud is 
optically thin to the ionizing radiation, as the models underpredict the amount of [Ne III] 
emission. Specifically, this result describes a matter-bounded component, in which some 
fraction of the ionizing radiation passes through the gas unabsorbed. Overall, a softer 
ionizing continuum will result in a lower [Ne V]/[Ne III] ratio, because there are fewer 
photons capable of ionizing Ne^^. 

Figure 6 shows that, independent of the ionization parameter, most of the observed 
ratios are matched by models with EUV/Soft X-ray spectral indices of 1.5 < ctj < 2.0. 
Using the 95% confidence interval for the actual mean for the [Ne V] 14.32^111 /[Ne III] 
and [Ne V]24.32^m/[0 IV] ratios observed in Seyfert 1 galaxies in our sample, we found 
a smaller range of spectral indices, i.e., 1.7 < cij < 2.0, with « 1.9 providing the best 
fit for the median observed values. A histogram of the EUV/Soft X-ray spectral indices, 

*The [Ne V]i4.32//m and [Ne III] fluxes were extracted from the Spitzer/IRS in the short- 
high order (SH, 4.7" x 11.3"). On the other hand, the [Ne V]24.32/im and [0 IV] emission 
were extracted from the Spitzer/IBS in the long-high order (LH, ll.l"x 22.3") 
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derived empirically from the position of each source in Figure 6, is presented in Figure 7. 
From their analysis of Hubble Space Telescope (HST) Faint Object Spectrograph (FOS) 
observations of a sample of 101 quasars at z > 0.33, Zheng et al. (1997) found that the 
power-law index in the EUV between 350 and 1050 A is o: = 1.96 ± 0.15. They also foimd 
that, for a sub-sample of 60 radio-loud sources, the EUV spectral index is a « 2.2, while 
for a sub-sample of 41 radio-quiet it is a « 1.8. As shown in Figure 6, most of the BAT 
AGN can be fitted with an EUV/Soft X-ray spectral index of a < 2.0, consistent with the 
radio-quiet nature of most Seyfert galaxies. 

From their analysis of HST spectra of a sample of 184 quasars at z > 0.33, Teller et al. 
(2002) found a slightly harder EUV continuum than Zheng et al. (1997). This sample 
included observations with the HST FOS, Goddard High Resolution Spectrograph eind 
available Space Telescope Imaging Spectrograph (STIS) data. For their complete sample, 
they found a spectral index of a = 1.76 ± 0.12 between 500 and 1200 A . For a subsample 
of radio-loud objects, they found oe = 1.96 ± 0.12, while for a radio-quiet subset, they foimd 
a = 1.57 ± 0.17. These values are roughly consistent with our range of EUV/Soft X-ray 
indices and the spectral index associated with the median values for the Se 3 dert 1 galaxies 
found in our sample, a as 1.9. Again, this supports the idea that a single power-law can 
represent the EUV/soft X-ray continuum in our sample of BAT AGN. One must note 
that star formation contamination in the [Ne III] emission can bias the EUV/Soft X-ray 
slope toweurds a steeper power-law index. In other words, star formation contamination 
in the [Ne III] emission will result in a lower [Ne V]i 4 . 32 /im/[Ne III] ratio, which could be 
wrongly associated with a softer EUV continuum. However, the tight correlation found in 
the [Ne III]-BAT relationship suggests that, on average, there is no strong enhancement 
due to star formation in the [Ne III] emission in our sample, although we cannot rule out 
relatively minor effects in individual objects. Despite the fact that our spectral index has 
been determined over a relatively small range of energy (~40-100 eV), these results support 
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the validity of our method, given the agreement with values derived from different spectral 
analyses, covering a wider range of energies. 

On the other hand, analysis of the Far Ultraviolet Spectroscopic Explorer {FUSE) 
archival data by Scott et al. (2004) revels a significantly harder composite spectral energy 
distribution, a = 0.56lo [If, than the spectral indices discussed above and no evidence of a 
UV spectral index break. One must note that, in their sample, some individual sources do 
show a spectral break in the far- ultraviolet (A ~ 910-2000 A). It should be noted that this 
is an heterogeneous sample of AGN, which were observed with FUSE because they were 
known to be bright in the UV. The authors reported a correlation between the spectral 
slope and AGN luminosity; in other words, low-luminosity AGN tend to show a harder 
ionizing continua. They suggested that the harder spectral index they derived is due to the 
relatively large fraction of low luminosity AGN in the FUSE sample®. In contrast, we found 
no evidence of such flatter indices in our sample. Specifically, in Figure 8 we compare our 
derived values for the spectral index and the 14-195 keV luminosities, from which it can 
be seen that there is no dependence between the EUV/Soft X-ray slope and luminosity in 
the BAT sample. Furthermore, Shang et al. (2005) found that most objects in their FUSE 
sample exhibit a spectral break around 1100 A , in agreement with' previous HST results 
(Zheng et al. 1997; Telfer et al. 2002). Their einalysis is based on a heterogeneous sample 
of 17 low-redshift AGN (z< 0.5), with quasi-simultaneous spectrophotometry including 
FUSE, HST and Kitt Peak National Observatory observations covering a rest wavelength 
from 900 — 9000 A. Given the fact that Scott et al. (2004) looked for possible systematic 
errors that could bias their results, the source of this discrepancy is imclear. While it is 
possible that the spectral break is luminosity dependent, another possible explanation is 

®The median value for their FC/5'E sample is log ALnoo = 45.0 (ergs/s) versus log ALnoo = 
45.9 (ergs/s) for the HST sample of active galaxies presented by Telfer et al. (2002) 
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that intrinsic reddening, if not fully corrected for, may lead to an incorrect determination 
of the spectral break and continuum shape (Shang et al. 2005). 

Figure 9 shows the predicted and observed [Ne V]i 4 . 32 /im/[Ne III] and [0 IV]/[Ne III] 
ratios. Contrary to the ratios presented in Figure 6, the [0 IV]/[Ne III] ratio could be 
susceptible to reddening and aperture effects. In this regard, these effects will have different 
manifestations in the observed [O IV]/ [Ne III] ratios. For example, the [O IV]/ [Ne III] 
ratio will have larger values if extinction is stronger for the shorter-wavelength line. On the 
other hand, aperture effects wiU tend to decrease the observed ratio. As seen in Figure 9, 
when using the [O IV] /[Ne III] ratio, we recovered the same range for the EUV/soft X-ray 
slope as those derived from the [Ne V]i 4 . 32 /xm/[Ne III] and [Ne V]24.32/im/[0 IV] ratios 
(see Figure 6). In other words, the [O IV]/[Ne III] ratios are fit by the range of spectral 
indices found from the other ratios, suggesting that reddening and aperture effects are not 
dominant (or, perhaps, cancel each other out) in the BAT sample of AGN, in agreement 
wth the discussion presented in Weaver et al. (2010). 


3.4. The AGN Parametrization 

Figures 10-11 show predictions for the [Ne V]i 4 . 32 /Ltm/[Ne III], [Ne V]24.32^m/[0 IV] 
and [0 rV]/[Ne III] ratios from a set of models generated with the parametrization of the 
AGN continuum including a UV bump, with cutoff temperatures from Tbb = 10®'^® K to 
10® ° K, and column densities of Nh = 10^^ cm~^ and Njj = cm“^. Once again, at 

Nh = lO^^cm”^ the cloud is optically thin to the ionizing radiation, thus underpredicts 
the [Ne III] emission. It is clear that the AGN parametrization with a UV “bump” 
temperature higher than 10®'^ K fails to predict most of the observed ratios, even in our 
radiation-bounded models. We find that the emission-line ratios for a majority of sources 
can be fit by models with UV bump temperatures of Tbb — 10®'^® K to 10®'^ K, or, 
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alternatively, peaks at ~ 600 1900 A. The best fit for the 95% confidence interval for 

the actual mean occurs between ~ 700 1000 A. These results are in agreement with the 

FUV break observed in the composite spectra of AGN, in the range of 1000-1100 A (e.g., 
Zheng et al. 1997; Telfer et al. 2002; Shang et al. 2005). 

This range of temperature is slightly lower than the typical temperatures used to 
characterize an average quasars continuum when modeling the BLR (Korista et al. 1997). 
However, the chosen form for the parametrization of the AGN continuum presented in 
Korista et al. (1997) focused on luminous quasars while our sample is dominated by lower 
luminosity sources. If the ionizing continuum is from an accretion disk, one would expect 
more luminous objects to be hotter for the same black hole mass, and, therefore, the lower 
range of temperature found in our sample is not unexpected. However, these results may 
also argue in favor of a different, intrinsic SED photoionizing the BLR. This intrinsic SED 
may be characterize by a bbb with a higher “blue bump” temperature, e.g., Tbb > 10®'^ K, 
than that found from our modeling of the NLR emission. In this interpretation, material 
located between the BLR and the NLR may absorb the intrinsic UV component. However, 
from our calculations, we do not find any clear evidence of a modification of the SED by an 
intervening medium. 


4. Conclusions 

In order to investigate the shape of the spectral energy distribution in the range 13.6 eV 
to 1 keV, we have generated photoionization models to characterize the small dispersion 
found in correlations between high-ionization mid-infraxed emission-line ratios in AGN 
This is the band where dust extinction and the line of sight absorber can make it virtually 
impossible to have a direct view of the continuum generated by the nuclear engine. Our set 
of radiation bounded models show that a broken power-law continuum model is sufficient 
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to reproduce the observed mid-infraxed [Ne V]i 4 . 32 /Ltm/[Ne III], [Ne V]24,32/im/[0 IV] 25 . 89 /im 
and [O rV] 25 . 89 /rm/[Ne III] ratios. Prom this we constrained the EUV/Soft X-ray spectral 
index, ctj, to be between 1.5 - 2.0, with ttj w 1.9 corresponding to the median value observed 
in Seyfert 1 galaxies. Our result agrees with previous studies (e.g., Zheng et al. 1997; 
Teller et al. 2002; Shang et al. 2005) and holds for the wide range of BAT luminosities and 
type I and II objects in the local universe, whereas previous studies of SED tend to focus 
on UV selected objects at a different redshift range. Despite the fact that our spectral 
index has been estimated for a relatively small range of energy (~40-100 eV), these results 
probe the validity of our method to uncover the AGN continuum in this energy bcind. Prom 
the comparison between the photoionization models and observations we derived empirical 
values for the EUV /Soft X-ray slope for each source and found no dependence between the 
EUV /Soft X-ray slope and the 14-195 keV luminosity in our sample. 

We also find that our power-law model can produce enough photons at energies >4 Ry 
to generate the amount of mid-infrared line emission observed in the BAT AGN sample. 
However, this result does not rule out the presence of a bbb. When we do add a “big 
bump” in the UV to an X-ray power-law, the majority of sources required “blue bump” 
temperatures of Tbb = 10®’^® K to 10®-^ K, suggesting that the peak of this component 
lies in a wide range of energies that can extend from ~ 600A to 1900 A. The best fit 
for the 95% confidence interval for the actual mean occurs between, ~ 700-1000 A, in 
agreement with previous estimates of the spectral index break (e.g., Zheng et al. 1997; 
Telfer et al. 2002; Shang et al. 2005). This range in temperatures is lower than the typical 
UV bump peak assumed to characterize an average quasar continuum when modeling the 
BLR (Korista et ad. 1997). However, if the ionizing continuum is from an accretion disk, 
one expects more luminous objects to be hotter for the same black hole mass. Therefore, 
the lower range of temperature found in our sample is not unexpected. 
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Fig. 1. — Comparison between the different continuum spectral energy distribution con- 
sidered in the present work. For the sake of comparison, for the “blue bmnp” model we 
present two values for the UV bump cutoff temperature 10® ™ K and 10® “ K. Again, for 
comparison purpose, we used a “power-law” model with a value for the EUV/Soft X-ray 
index of Oj = 1.5. The three continua shown are normalized to have the same integrated 
luminosity between 1 Ryd< hu < 7.354 x 10® Ryd, however, because we are focusing on the 
modeling of emission line ratios, instead of of absolute line fluxes, the chosen normalization 
is arbitrary. 
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Fig. 2. — Upper panel. The [Ne V] ratio versus the hydrogen number density, Uh for 
the “power-law” model with ccj = 1.5 for different ionization parameters, logU = —0.5 
and log U = —3.0. Lower panel. The [Ne V] ratio versus the hydrogen number density. 
Tiff, for the “power-law” model for different spectral slopes with a fix ionization parameter, 
logU = —1.5. The solid horizontal line represents the observed average ratio and the 1-(t 
{dashed lines) found for Seyfert 1 galaxies. 
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Fig. 3. — Predicted ionic column densities for Ne III (circles), Ne V (triangles) and O IV 
(squares), for two set of models, as a function of ionization parameter U, and two values of 
the column density, lO^^cm"^ and lO^^cra"^, left and right panels, respectively. The solid 
line represents the “power- law” model with Oj = 1.5 and the dashed line represents the “big 
burpp” parametrization with the peak of the UV bump at w 44 eV, or alternatively a UV 
bump cutoff temperature of 10® K. 
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Fig. 4. — Comparison for the radiation bounded models between the total Hydrogen col- 
umn density and ionization parameter for the “power-law” model for different values of the 
EUV/Soft X-rav index, cci. One must note from this comparison that, within the range 
of ionization parameters that reproduce the observed ratios, our radiation bounded mod- 
els yield to Hydrogen column densities higher than that found in typical NLR conditions, 
logNff = 21 {cm“^) (e.g., Kraemer et al. 2000) 
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Fig. 5. — A comparison between the ionization parameter associated with the peak of the 
ionic column density versus the EUV/Soft X-ray spectral slopes for the “power-law” model. 
Prom this comparison is clear that the peak on the ionic colunm density depends on the 
shape of the ionizing continuum. 
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Fig. 6. — Comparison betv.'een the predicted and observed [Ne Vji^ 32 jtitn/[Ne III] and 
[Ne V]24.32/im/[0 IV] ratios. Seyfert 1 galaxies are presented as black circles, Seyfert 2 galax- 
ies are red triangles, blue squares represent the newly detected BAT AGN and green stars are 
LINERs. We present two set of “power-law” models corresponding to logiVA = 21 (cm“^) 
{solid lines) and logNn = 22 (cm~^) {dashed lines) for a range of ai — 1.5 — 2.0. For the 
set of models presented in this comparison we show a range of ionization parameters of 
—2.5 < log 17 < —1.5, in steps of 0.5 dex and 0.1 dex for logiV^ = 21 and logAT^ = 22, 
respectively. All these models corresponds to a hydrogen number density of = 10^ (cm®). 
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Fig. 7. — Distribution of EUV/Soft X-ray slopes, ai, derived from our analysis of the entire 
BAT sample. One must note that the spectral index was derived empirically from the 
position of each source in Figure 6, therefore, each slope has an error. For any given 
source located between two consecutive power-law models, e.g., different spectral slope, we 
selected the flatter index as the empirical estimate. 
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Fig. 8. — Comparison between the spectral index, ai, and the 14-195 keV luminosities 
where it can be seen that there is no dependence between the EUV/Soft X-ray slope and the 
14-195 keV luminosity in our sample. One must note that the spectral index was derived 
empirically from the position of each source in Figure 6, therefore, each slope has an q,±1 
error. For any given source located between two consecutive power- law models, e.g., different 
spectral slope, we selected the flatter index as the empirical estimate. 
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Fig. 9. — Comparison between the predicted and observed [Ne V]x 4 . 32 /zm/[Ne III] and 
[0 rV]/[Ne III] ratios. Seyfert 1 galaxies are presented as black circles, Seyfert 2 galaxies 
are red triangles, blue squares represent the newly detected BAT AGN and green stars are 
LINERs. We present two set of “power-law” models corresponding to log Njj = 21 (cm“^) 
{solid lines) and logiV/x = 22 (cm“^) {dashed lines) for a range of Oi = 1.5 — 2.5. For the 
set of models presented in this comparison we show a range of ionization parameters of 
—2.5 < log 17 < —1.0, in steps of 0.5 dex and 0.1 dex for logNu = 21 and log Njj = 22, 
respectively. AH these models corresponds to a hydrogen number density of Uh = 10^ (cm^). 


30 



Fig. 10. — Comparison between the predicted and observed [Ne V]i 4 . 32 ;um/[Ne III] and 
[Ne V]24.32//m/[0 IV] ratios. Seyfert 1 galaxies are presented as black circles, Seyfert 2 
galaxies are red triangles, blue squares represent the newly detected BAT AGN and green 
stars are LINERs. We present two set of “blue bump” models corresponding to log Ng = 
21 (cm“^) (solid lines) and log A// = 22 (cm“^) (doshed lines) for five values of the UV bump 
cutoff temperature between Tbb = 10®‘^® K and 10® ** K. For the set of models presented 
in this comparison we show a range of ionization parameters of —2.5 < logt/ < —1.0, in 
steps of 0.5 dex for both logNff = 21 and log Njj = 22. All these models corresponds to a 
hydrogen number density of = 10® (cm®). 
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Fig. 11. — Compeirison between the predicted and observed [Ne V]i 4 . 32 /xm/[Ne III] and 
[C» IV] /[Ne III] ratios. Seyfert 1 galaxies axe presented as black circles, Seyfert 2 galaxies 
are red triangles, blue squeures represent the newly detected BAT AGN and green stars are 
LINERs. For the sake of simplicity, we present two set of “blue bump” models corresponding 
to logNji = 21 (cm“^) {solid lines) and logNjj = 22 (cm~^) {dashed lines) only for three 
values of the UV bump cutoff temperature, Tbb — 10®'^® K, 10® ™ K and 10® ® K. For the 
set of models presented in this comparison we show a range of ionization parameters of 
—2.5 < logC/ < —1.0, in steps of 0.5 dex for both logiV^ = 21 and logAT^ = 22. All these 
models corresponds to a hydrogen number density of n* = 10® (cm®). 


Table 1. 


Statistical Analysis Between Seyfert 1 and Seyfert 2 Galaxies 



Seyfert 1 


Seyfert 2 


Pk-3 

(%) 

Measurements 

Available 

Mean 

Standard 

Deviation 

Measurements 

Available 

Mean 

Standard 

Deviation 

[Ne V]i 4 . 32 ^rn/[Ne V]i4.32,im 

29 

0.95 

0.28 

30 

1.00 

0.38 

87.8 

[Ne V]i 4 . 32 #*m/[Ne HI] 

32 

0.47 

0.22 

33 

0.41 

0.20 

79.8 

[Ne V):4.3j/xni/[0 IV| 

31 

0.34 

0.19 

30 

0.28 

0,06 

57.0 

(O rv]/[N6 III] 

38 

1.58 

0.87 

33 

1.65 

0.77 

91.4 


Note. — The last column, Pk-s» represents the Kolmogorov-Srnirnov (K-S) test null probability. Upper limits for the 
[Ne V] fluxes are not included. This table also includes information about the numbers of Seyfert 1 and Seyfert 2 galaxies, 
mean '*a.ues and standard deviations of the mean for the measured quantities. For the [Ne V]i 4 . 32 /im/[Ne HI] ratios the 
95% confidence interval for the actual mean is 0.3904 - 0.5465 for Seyfert 1 galaxies, and 0.3417 thru 0.4831 for Seyfert 2’s. 
For the [Ne V]24.32/xm/[0 IV] ratios the 95% confidence interval for the actual mean is 0.2717 thru 0.4089 for Seyfert 1 
gala:des and 0.2469 thru 0.2964 for Seyfert 2’s. Finally, for the [O IV] /[Ne III] ratios the 95% confidence interval for the 
actual mean is 1.295 thru 1.867 for Seyfert 1 galaxies, and 1.375 thru 1.924 or Seyfert 2’s. 
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